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ABSTRACT 

GRB 130925A is an ultra-long GRB, and it shows clear evidences for a thermal emission in the soft X- 
ray data of Swift/XKT 0.5 keV), lasting till the X-ray afterglow phase. Due to the long duration of 
the GRB, the burst could be studied in hard X-rays with high-resolution focusing detectors (NuSTAR). 

The blackbody temperature, as measured by the Swift/XRT, shows a decreasing trend till the late 
phase (Piro et al. 2014) whereas the high-energy data reveals a significant blackbody component 
during the late epochs at an order of magnitude higher temperature (~ 5keV), as compared to the 
contemporaneous low energy data (Bellm et al. 2014). We resolve this apparent contradiction by 
demonstrating that a model with two black bodies and a power-law (2BBPL) is consistent with the 
data right from the late prompt emission to the afterglow phase. Both the blackbodies show a similar 
cooling behaviour upto the late time. We invoke a structured jet, having a fast spine and a slower 
sheath layer, to identify the location of these blackbodies. Independent of the physical interpretation, 
we propose that the 2BBPL model is a generic feature of the prompt emission of all long GRBs, and 
the thermal emission found in the afterglow phase of different GRBs reflects the lingering thermal 
component of the prompt emission with diverse time-scales. We strengthen this proposal by pointing 
out a close similarity between the spectral evolutions of this GRB and GRB 090618, a source with 
significant wide band data during the early afterglow phase. 

Subject headings: gamma-ray burst: general — gamma-ray burst: individual (130925A) — methods: 

data analysis — methods: observational — radiation mechanisms: thermal 


1. INTRODUCTION 

The X-ray and gamma-ray emissions from gamma-ray 
bursts (GRBs) show two distinct phases: a prompt emis¬ 
sion with drastic spectral and temporal variations, and 
an aft erglow havin g smooth variation with time and en¬ 
ergy (jPiranI Il999f) . Though the afterglow emission is 
understood as due to synchrotron emission from a fast 
moving jet in the ambient medium, the characteristics 
and the radiation mecha nisms of the pro mpt emission is 
still a matter of debate (|Meszarosll2006ll . This apparent 
dichotomy of observations, and the corresponding un¬ 
derstanding, is accentuated by the fact that the prompt 
emission is studied using wide band hard X-ray/gamma- 
ray detectors with poor energy resolution and low sensi¬ 
tivity (due to high background), while the leisurely X-ray 
afterglow is amenable to be observed with focusing X-ray 
telescopes with excellent energy res olution and sensitiv- 
ity. The fast moving Swift satellite (jGehrels et al.l[20Ml 
provides good quality data for the early afterglow of sev¬ 
eral G RBs in the low energ ies with the X-Ray Telescope 
IXRT: [Burrows et al.l 1200511 . However, the typical slew¬ 
ing ability (~ 60 — 80s) is insufficient to capture the 
prompt emission of the majority of GRBs with the XRT. 
Further, during the XRT observations, the hard X-ray 
emission fades below t he detection level of Sw ift/Burst 
Alert Telescope (BAT: [Barthelmv et al.l[2005fl for most 
cases. Hence, it is difficult to identify a broadband spec¬ 
tral feature with th e combined BAT-XR T data. 

In a recent paper (jBasak fc R,aol[^I4aL hereafter Paper 
I), we found a clear evidence for two smoothly evolving 
blackbodies in the early afterglow of GRB 090618. As 


this GRB has a long gamma-ray duration, a significant 
overlap is seen in the BAT and XRT observations dur¬ 
ing the junction of the prompt and the afterglow phase 
(125 — 165 s). This overlapping observation allowed us 
to identify two blackbodies in the individual energy win¬ 
dows provided by the BAT and the XRT. The temper¬ 
ature and the flux of the two blackbodies are found to 
have similar evolution throughout the final pulse of the 
prompt emission till the early afterglow phase. This is 
one of the rare GRBs with a significant overlap in the 
XRT and BAT observations, and further, the two black¬ 
bodies happen to independently dominate the two detec¬ 
tors. Hence it is important to confirm this finding and 
examine whether this behaviour is a generic feature of 
all GRBs. 

In this context ultra-long GRBs, a class of GRBs which 
produce prompt emission for hours, provide an excel¬ 
lent opportunity to study the spectral evolution from 
the prompt to the afterglow phase. GRB 130925A is 
an ultra-long GRB with the gamma-ray duration of 
^ 2 hr. This GRB was detected by several satellites in¬ 
cluding the Swift/BAT and the Swift/ABT. Due to the 
gamma-ray emission on a long time-scale, the burst could 
be observed with th e focusing telescopes of NuSTAR 
((Harrison et al.l[2M^ . The Swift/ART data revealed the 
presence of a blackbody emission a t low energies du ring 
the prompt to the afterglow phase ()Piro et al.ll2ni4 P14 
hereafter), though the same data could be interpreted as 
due t o a steep power-la w emission rather than a black¬ 
body ((Evans et al.ll20l3l . The hard X-ray data from the 
NuSTAR, on the other hand, found a very significant 
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emission above ~10 keV at late times, which could be 
interpreted as due to an emission from a blackbody of 
temperature ~5 keV (|Bellm et alJl20ll: B14 hereafter). 
The contemporaneous data of Swift/XRT, however, finds 
a blackbody with temperature ^0.5 keV (P14). In this 
Letter we demonstrate that all these data are consistent 
with a model consisting of two black bodies and a power- 
law (2BBPL). In the next section we summarise the per¬ 
tinent data on GRB I30925A and present a re-analysis. 
The results are discussed in the context of a spine-sheath 
jet (section 3), and in the last section we give a summary 
of the important hndings of this work. 

2. THERMAL EMISSION IN GRB 130925A 

GRB 130925A was triggered by the Swift/BAT at Tq = 
2013-09-25 04:11:24 UT. The Swift/XRT followed the 
burst from 150 s after the BAT trigger. The burst showed 
several pulse emissions and soft flares till late time. A 
detailed analysis of th e high-energ y emis sion from this 
source is presented in lEvans et al.l (|2014[1 . They argue 
that this GRB (and by extension other ultra-long GRBs) 
could be understood as events with low circumburst den¬ 
sities, and correspondingly having an order of magnitude 
lower deceleration time compared to the normal long 
GRBs. They strongly favour a steep power-law function 
(index P ^ —4) for the late X-ray afterglow spectrum, 
and use the dust scattering model for the explanation 
(IShao fc DaillMlTll . 

P14, on the other hand, explain the steep spectrum 
as a superposition of a blackbody (temperature, kT ~ 
0.5 keV) and a power-law (BBPL) model. A crucial 
XMM-Newton observation after 3 months shows a harder 
power-law index (~ —2.5). As the thermal emission pos¬ 
sibly subsides after such a long time, the spectrum is 
expected to reveal the underlying non-thermal compo¬ 
nent. The hard value of the index agrees with a stan¬ 
dard afterglow in a wind medium. On the other hand, 
the blackbody component is suggested to be a ‘hot’ co¬ 
coon which forms on the GRB jet as it pierces through 
the progenitor — possibly a blue suergiant star. 

Another important clue on the spectral shape is pro¬ 
vided by the observation with high-resolution focusing 
detectors, NuSTAR, in three periods (1.8 day, 8.8 day, 
and 11.3day), and Chandra AGIS-S in one period 
(11.0 day). B14 have analyzed the spectra in two epochs 
— ^ 2 day and ^ 9 — 11 day. In both cases, they 
have found significant deviation from the canonical af¬ 
terglow spectrum. The high-resolution data shows a 
‘dip’ in the continuum spectrum. B14 have used a 
BBPL model to capture the feature. Interestingly, they 
have found that the blackbody temperatures in the two 
epochs (5.6lj;2keV and 4.0lQ;gkeV, respectively) have 
an order of magnitude higher values than that found 
by P14 in the XRT data. Though the data allows a 
softer blackbody (near ^ 0.5keV), B14 have argued in 
favour of the higher-temperature blackbody, as the lower- 
temperature blackbody shows an unphysical contraction 
of the apparent emitting region at the later epoch, from 
(3.2 ± 0.8) X 10^° cm to (1.7 ± 0.4) x 10^° cm. 

2.1. Signature of two blackbodies in the afterglow data 

To arrive at a correct spectral model, it is important 
to hx the continuum and hence data above ^10 keV is 
very crucial. We re-analyze the data whenever a higher 


Table 1 

Parameters of the two blackbodies in the two NuSTAR 
observation 


Interval 

kT 

Pbb 

Rbb 


(keV) 

(ergcm“^ s“^) 

10° cm 

Observation 1 Tq + 1.8 day) 

G)BBi 


(5.46 ±2.58) X 10-1° 

0.9 ±0.3 

BB2 

0.421°;“ 

(1.25 ± 0.92) X 10-“ 

(2.0 ± 1.0) X 10° 

Observation 2 To + 11 day) 

BBi 

o q+0.6 

(1.3 ± 0.3) X 10-1° 

1.3 ±0.2 

BB2 

0.15 ±0.03 

(2.0 ± 1.4) X 10-“ 

(2.0 ± 1.0) X 10° 


(‘^)BBi=Higher-temperature blackbody, BB 2 =Lower-temperature 
blackbody 


energy observation is available (either BAT or NuSTAR 
observation). These epochs are: (I) 150 — 300s dur¬ 
ing the prompt emission (joint BAT-XRT observations), 
(II) 1.8 day during the afterglow (joint NuSTAR-XTTT 
observations), and (III) ~ 9 — 11 day (joint Chandra- 
NuSTAR observations). The observation in epoch II 
provides the best possible data among all the observa¬ 
tions, due to the superior sensitivity and resolution of 
the NuSTAR detectors (and the higher flux compared 
to epoch III). In Figure [U we show the spectrum of 
epoch II, fitted with various models. We follow the stan¬ 
dard procedure described by B14 to extract the spectral 
data. The spectrum is fitted in XSPEC vl2.8.2. We fit 
the following models: (A) power-law (PL), (B) Power- 
law with a Gaussian absorption (gabsxPL), (G) a black¬ 
body with a power-law (BBPL), and (D) two blackbod¬ 
ies with a power-law (2BBPL). In all cases we have used 
a wabs and a zwabs model to account for the absorp¬ 
tion in our Galaxy and in the source, respectively. P14 
have fitted the time-resolved data during 20 ks-8 Ms with 
wabs X zwabs X PL model. As the equivalent hydrogen col¬ 
umn density (Vh) at the source should not change over 
time, they have linked this parameter in all time bins to 
obtain its precise value — (2.1 ± 0.1) x 10^^cm“^. We 
have used this value, along with the Galactic absorption 
of 1.7 X 10^° cm“^ to perform uniform fitting with all the 
models. Keeping the source absorption free gives a value 
close to P14, within the statistical errors. Note in Fig- 
ure[l] (panel A) that a PL fit shows a ‘dip’ near 5 keV. We 
obtain degrees of freedom (d of)=158.7/124 for a PL 
fit. Incidentally, a Band function (iBand et al.lll993ll gives 
a similar unacceptable fit with x^/dof=153.0/121. A 
Gaussian absorber gives an immediate improvement with 
y^/dof=100.8/121. However, as noted by B14, we also 
see a shift in the centroid from 6.4±0.8keV to < 4.0 keV 
from epoch II to epoch III. Hence, the Gaussian absorber 
is unphysical. A BBPL model gives a reasonable fit with 
X^/dof=106.1/122 (see panel G). Finally, we fit the spec¬ 
trum with a 2BBPL model and obtain x^/dof=96.7/120. 
The parameters of the blackbodies are shown in the first 
two rows of Table [TJ We perform F-test to find the sig¬ 
nificance of adding the second blackbody. We find a 2.9cr 
significance corresponding to p —value = 3.8x10“^. The 
temperature of the two blackbodies in this epoch are: 
5.12l^'52keV and 0.42^°keV. Note that these values 
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Table 2 

Parameters of the BBPL model fitted to the BAT data in 105-150 s interval, and that of the 2BBPL model fitted to the joint BAT-XRT 

data in 150-295 s interval 


Parameter 

105-150s 

150-175 s 

175-225s 

225-295 s 

riH ( 10 ^^ atoms cm“^) 


2.1 ± 0.2 

1.7 ±0.2 

1.7 ±0.2 

kT (keV) 

13.0/^;° 




N 

nil +0.09 




(keV) 


8.91/® 


S-Stru 

Nh 


10 

0 q 

d d 

+ 1 

cd 

(3.5 ±3.2) X 10-2 

(3.2t|j) X 10-2 

kTi (keV) 


9 t+0.5 
^•-*-- 0.6 

2.1 ±0.5 

9 q + 0.5 
^•'^-0.7 

Ni 


0 14+0-07 

^•-L^_0.06 

X 10-2 

X 10-2 

r 

- 2.3 ± 0.2 

- 2.2 ± 0.1 

- 2.1 ± 0.1 

- 2.1 ± 0.1 

X?ed (dof) 

0.93 (54) 

1.2 (76) 

1.3 (76) 

1.1 (76) 


('^) /i=Higher-temperature blackbody, Z=Lower-temperature blackbody 



Energy (keV) 


Figure 1. Spectral fitting to the combined data of NuSTAR and 
S'm/t/XRT of GRB 130925A during 1.8 day after the S'm/t/BAT 
trigger (see text for explanation). 

are close to those reported by B14 and P14 (both within 
Icr), respectively. Based on the analysis of the epoch II, 
we fit the spectrum of the epoch III with a 2BBPL model. 
We freeze the power-law index to the value obtained in 
epoch II (-3.74). B14 has combined the NuSTAR data of 
To + 8.8 day and Tq -I-11.3 day. As both the blackbodies 
may evolve significantly during this time, we use only the 
data of To -I- 11.3 day. We obtain a x^/dof=146.55/133. 
The parameters of the blackbodies are shown in the last 
two rows of Table [H The temperature of the two black¬ 
bodies are 3.0to5keV and 0.15±0.03keV. We note that 
the temperature of the higher-temperature blackbody is 
somewhat lower than the temperature of the single black¬ 
body as obtained by B14 (d.Ot^o). 

2.2. Thermal emission spanning the prompt and 
afterglow phase 

In order to see how the temperature of the two black¬ 
bodies evolve during the prompt emission, we fit the joint 
BAT-XRT data in 150-300 s interval. We obtain three 
time-resolved spectra in 150-175 s, 175-225 s, and 225- 
295 s intervals. For each case, we obtain a reasonable fit 
with reduced close to 1 for 76 dof (see Table [21). In 
the first interval, the 2BBPL model shows significant im¬ 


provement compared to a BBPL model with = 18.6 
in the expense of two dof. The F-test gives a 3.3cr sig¬ 
nificance {p — value = 9.3 x 10“"'^) for the addition of the 
second blackbody. It is worthwhile to mention that a 
Band function also gives an acceptable fit to the prompt 
emission data. However, as already mentioned the Band 
function is inadequate during the epoch II, where we have 
good quality data. In addition to the joint BAT-XRT 
data, we also ht a BBPL model to the falling part of the 
final pulse (105-150 s) and obtain kT = I 2 . 91/3 Y keV. 

In Figure m we show the temperature evolution of the 
two blackbodies (filled symbols — circles for the higher-, 
and squares for the lower-temperature blackbody). The 
lightcurve (15-50 keV flux in units of ergcm“^ s“^) of the 
GRB in the BAT and XRT detectors is shown in the 
background. The star represents the single blackbody 
temperature during the falling part of the last pulse in 
the BAT data. For a comparison, we have plotted the 
temperature values obtained by P14 (open squares) and 
B14 (open triangles). We note that the values of the 
higher and the lower temperature are in general agree¬ 
ment with that of B14 and P14, respectively. We fit 
the evolution of both the blackbodies as a power-law 
function of time. In doing this, we assume that all 
the blackbody temperatures given by P14 are that of 
the lower-temperature blackbody of the unified 2BBPL 
model. The slope of the evolution are —0.15 ± 0.05 and 
—0.20±0.02 for the higher- and lower-temperature black¬ 
body, respectively. The two evolutions can be considered 
either approximately similar, or the higher-temperature 
blackbody has a slower temperature evolution. Note that 
due to the large gap in the prompt and afterglow cov¬ 
erage of the high energy observation, the error in the 
slope of the higher-temperature blackbody evolution is 
large. Hence, the comparison between the two temper¬ 
ature evolutions remains inconclusive. In the inset of 
Figure m we show the time evolution of the flux of the 
individual components of the 2BBPL model. We have 
shown only those points, where we have high energy ob¬ 
servation with either BAT or NuSTAR. In order to com¬ 
pare the flux evolution of the components, the blackbody 
flux are normalized to the power-law flux of the initial 
time bin. We note that the thermal flux decreases more 
rapidly compared to the power-law flirx. The flux of the 
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Figure 2. Evolution of the temperature of the two blackbod- 
ies found in various epochs of GRB 130925A (filled symbols for 
present work, star for the falling part of the last pulse — 105- 
150 s). The data from P14 and B14 are shown by open squares, 
and open triangles, respectively. The 15-50 keV flux is shown in 
the background. Inset: Flux evolution of different components — 
filled circles: higher-temperature blackbody, open circles: lower- 
temperature blackbody, plus: power-law. The flux of the two 
blackbodies are normalized with the power-law flux at the initial 
time bin (multiplied by 13.6 and 10.7, respectively). 


higher-temperature blackbody evolves even faster. 

3. DISCUSSION 

Though it is not conclusive that the 2BBPL is the cor¬ 
rect description of the GRB emission model during the 
prompt and the early afterglow phase, there are strong 
evidences to support the claim that this model is consis¬ 
tent with a wide variety of data. This model is superior to 
other models (BBPL/ Band) for bright GRBs (|Rao et al.l 
I2014D and it is a preferred mo del for GRBs with sin¬ 
gle pul ses (iB asak fc R,aoll2pi4bt) . GRBs with separable 
pulses (jBasak fc R,aoll2013b[ ) and GRBs with high energy 
(GeV) emission and detecte d by Fermi/Large Area Tele¬ 
scope (|Basak fc Raoll2013^ . For GRB 090618 (Paper I), 
the superior spectroscopic Swift/XRT detector shows a 
blackbody temperature decreasing with time, consistent 
with the interpretation that the two components of the 
2BBPL model seen in the last GRB pulse are showing 
a cooling behaviour. In the present paper, we have seen 
yet another evidence for the two blackbody emissions in 
the case of GRB 130925A. In the next sections, we in¬ 
vestigate the implications for the putative photospheres 
and a possible jet scenario where two photospheres can 
co-exist. 


3.1. Photosphere of the two blackbodies 

It is interesting to calculate the apparent emission ra¬ 
dius for a spherical emission region (Abb) during the 
two epochs of NuSTAR observation. Both 1314 and P14 
find that the value of i?BB shows a contraction during 
this time. As we have two blackbodies, we calculate the 
value of Rbb for each of them as follows. 


\/(t[(1 + z)T]4 

where dg is the luminosity distance, a is Stefan- 
Boltzmann constant = 5.6704 x 10“® ergcm“^ s“^ K“'*, 
Fbb is the observed blackbody flux, and (1 + z)T is 


the temperature corrected for the redshift, z. The ra¬ 
dius of the higher-temperature blackbody in the two 
epochs are found to be Rbb = (0.9 ± 0.3) x 10 ® cm and 
Rbb = (1.3±0.2) x 10 ® cm, which show a slow increment. 
Interestingly, the value of Abb of the lower-temperature 
blackbody is found to be Rbb = ( 2.0 ± 1 . 0 ) x 10 ^° cm 
and i?BB = ( 2.0 ± 1 . 0 ) x 10 ^^ cm (see the last column 
of Table [T|). Note that this result is markedly different 
from that obtained by B14 and P14. Our analysis shows 
a definite increment of the apparent radius at the later 
epoch. 


3.2. A spine-sheath jet 

In paper I, we have proposed a spine-sheath jet to ex¬ 
plain the evolution of the two blackbodies. A fast moving 
spine with a slower sheath component is expected in a 
number of physical scenario, e.g., a hot cocoon formed 
over the G RB jet as it pierces through the envelop of the 
progenitor (iMeszaros fc R.eesI 2001b iR.amirez-R.iiiz et al.1 
l200a IZha.ns- et a.l.l l2003L I2004l)~ or a collimated proton 
j et along with a wider neutron sheath fo r a magnetic jet 
(|Vlahakis et al.l[^003b iPeng et al.l 120051) . Such a struc¬ 
tured jet is frequently invoked to explain double jet 
break and optical re-brightening (e.g..lBerger et al.l[2003b 
iLiang fc Pail l2004b iHolland et al.l 120121) . If the sheath 
component of the jet is in fact the cocoon, we expect 
it to be less collimated, and the terminal bulk Lorentz 
factor ( 77 ) to be lower by a factor of 10 than the spine. 

The observed physical sizes of the two blackbodies 
show very slow expansion rates and the correspond¬ 
ing sizes of the photospheres too would have slow 
speeds and hence they cannot be tied to the actual 
sizes of the outflowing plasma. Hence the photospheres 
have to be steady/ gradually changing regions in the 
outflow. We invoke that the spine photosphere pro¬ 
duces the higher-temperature blackbody while the lower- 
temperature blackbody is produced by the sheath photo¬ 
sphere. Note that the size of the emitting regions calcu¬ 
lated above are considered as spherical emission regions. 
The physical photosphere can be related to the observed 
size by rph = pRbb- With reasonable values of rj for the 
spine and the sheath (say 10 ® and 10 ^, respectively), one 
should obtain larger values of the corresponding emis¬ 
sion radius. The observed properties of the two photo¬ 
spheres can be used to constrain the properties of the 
spine and the sheath. For example, we have found that 
the higher-temperature blackbody has much lower evo¬ 
lution and hence we can assume that the spine remains 
“steady” for a long time. The sheath photosphere, on 
the other hand, shows considerable evolution. For a co¬ 
coon shea th, a rapid photosp heric expansion is indeed 
expected (|Starling et al.l[2M^ . 

3.3. Comparison with GRB 090618 

lEvans et al.l (I2014D have made a systematic analysis of 
672 XRT GRBs, and found that after 3 ks of the indi¬ 
vidual trigger, the hardness ratio (HR) does not show 
any evolution for the majority of GRBs. For 12 GRBs 
(including GRB 130925A), which show > 5cr significance 
of HR variation at late times, the power-law index shows 
a continuous increment with time. They argue that the 
HR evolutions in these GRBs including GRB 130925A 
can be explained by a dust scattering model. 
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Figure 3. A comparison of spectral evolution between the ultra- 
long GRB 130925A with GRB 090618, a typical long GRB (open 
and filled symbols, respectively). The 100-300s time axis of the 
latter is stretched in the interval 100 s to 2 X 10^ s. Upper panel: 
Temperature evolution of the two blackbodies. We have shown the 
15-50 keV flux in the background (higher flux for GRB 090618). 
Lower panel: Hardness ratio (count ratio in 1.5-10 keV to that of 
0.3-1.5keV) of the two GRBs. 

However, the hard X-ray data at late times requires 
emission above 10 keV, which is not compatible with a 
steep power-law model (as required by the XRT data 
alone). Hence, we explore here whether the spectral evo¬ 
lution of GRB 130925A could be explained by the 2BBPL 
model, and whether this model is a generic feature of 
all GRBs. We have taken the temperature evolution as 
given in Figure [21 and verified that the count rate varia¬ 
tion in two energy bands in the XRT data (0.3 — 1.5 keV 
and 1.5 — lOkeV) after 100 ks is compatible with the flux 
variation of the three components as given in the inset 
of Figure m 

The rare occurrence of spectral variation after 3 ks 
could be explained by postulating that the blackbody 
components decay at different rates for different GRBs. 
We investigate this hypothesis by comparing the results 
of GRB 130925A with GRB 090618. In Figure|31 we plot 
the relevant parameters for these two GRBs, but the time 
axis of GRB 090618 is arbitrarily stretched (100 — 300 s 
data is stretched to 100 — 2xl0^s). The close similarity 
in the spectral variation in these two GRBs encourages 
us to hypothesise that the 2BBPL model is a generic fea¬ 
ture of all GRBs and the black body components decay 
at different rates for different GRBs. 

4. CONCLUSIONS 

The major conclusions from the work presented here 
can be summarized as follows. 

(1) We use high-resolution NuSTAR data of the ultra- 
long GRB 130925A, and find that the spectrum is consis¬ 
tent with a 2BBPL model. The lower-temperature black- 
body is consistent with an evolving blackbody through¬ 
out the burst (P14). The higher-temperature blackbody 
is also consistent with the single blackbody of B14, and 


this component also shows a similar evolution as the 
lower-temperature blackbody. 

(2) The evolution of the corresponding photospheres of 
the two blackbodies are found to be different. While the 
higher-temperature blackbody shows a constant spher¬ 
ical emission region, the lower-temperature blackbody 
shows a rapid evolution of the photosphere. 

(3) We have proposed a structured jet with fast spine 
and slow sheath as the origin of the thermal components. 
The higher-temperature blackbody is produced at the 
spine photosphere, which remains ‘steady’ at the later 
epoch. The sheath is possibly a cocoon, which shows a 
rapid expansion of the photosphere. 

In recent years, we have found that the time-resolved 
prompt emission spectrum of a variety of GRBs is co n- 
sistent with the 2BBPL model (iBasak fc Raol I^13bllal : 
IRao et aDl2014D . This spectral shape is statistically pre¬ 
ferred for GRBs with high signal-to-noise data. However, 
due to the limited spectral resolution of the gamma-ray 
detectors, we could not draw a firm conclusion. We could 
detect the two blackbodies in the early afterglow data 
of GRB 090618 (Paper I), as the XRT provided a good 
spectral data during the overlapping observation with the 
BAT. The current observation also shows a clear signa¬ 
ture of the two blackbodies which evolve from the very 
early phase during the prompt emission, and remains vis¬ 
ible with high significance at very late times. We empha¬ 
size again that independent of the physical interpretation 
of the 2BBPL, this model emerges as a generic spectrum 
of long GRBs. 
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